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Abstract: This is an expansion of a talk for mathematics and physics students of the 
O , Manchester Grammar and Manchester Hisrh Schools. It deals with numbers such as 

the Riemann zeta value ({3) = J2n>o ^/n ■ Zeta values appear in the description of 
P^l sunshine and of relics from the Big Bang. They also result from Feynman diagrams, 

■ which occur in the quantum field theory of fundamental particles such as photons, 

■ electrons and positrons. My talk included 7 reasonably simple problems, for which 
^ . I here add solutions, with further details of their context. 
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1 Numbers 

1.1 Sum of the first 136 cubes 



o 
o 

1^ = 1 = l2 

1^ + 2^ = 9 = (1 + 2)2 

13 + 2^ + 3^ = 36 = (1 + 2 + 3)2 

>: 1^ + 2^ + 3^ + ... + 136^ = (68 X 137)2 = (1 + 2 + 3 + ... + 136)2 



It is rather easy to prove, by induction, that 



OO 

cn 

^ ■ N-l . 

O ■ n=l 



for every integer N > 2. 

Assume that formula ([1]) is true for one particular integer, say N = M. Then 



U- V = -M2(M - 1)2 + M' = -M\M + 1)2 

„=i 4 4 

and ([1]) is true for A = M + 1. It is true, by inspection, for A = 2, and hence for 
every integer A > 1. 

1.1.1 Problem 1: proof by induction 

Prove that, for every integer A > 1, 

N-l 1 

^ = _Ar(Ar- 1)(2A- 1)(3A2 -3A- 1). (2) 

n=l 30 
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1.2 Sum of the first 136 inverse cubes 



n=l 



E 




3^ 
23 



n=l 




23 X 33 



251 




5 X 11 X 37 
26 X 33 

72 X 11 X 53 



29 X 3 X 53 X 73 



29 X 36 X 53 X 73 X 113 



11^ X 89 X 359 X 4957 
29 X 36 X 53 X 73 

13^ X 151099201553 



26 X 3 X 53 



31 X 2538983 



It is quite hard to obtain the complete factorization of the numerator of the sum of 
the first 136 inverse cubes. Let 



then the numerator of 5" (137) turns out to be 

137^ X 359 X 20273371 x 20077753681 x 95441340948666564707 x P29 x P47 x P54 
where 

P29 = 11670683543184939914296499797 

P47 = 15515586948180802047607109239107700793654018233 

P54 = 135756874001680563462725618538328344771913356122240503 

are primes with 29, 47 and 54 decimal digits. 
1.2.1 Problem 2: modular arithmetic 

Find all the primes p such that divides the numerator of S{p), using the result 
that for every prime p and every integer n with p > n > there is a unique integer 
m with p > m > and mn — 1 divisible by p. [Hints follow.] 



N-l 



1 

n3 



S{N) = E 



n=l 
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1.2.2 Hints for Problem 2 



Let p be any prime. Let a and b be rational numbers with denominators that are 
not divisible by p. If the numerator of a — 6 is divisible by p, we say that a and b 
are congruent modulo p. Let's agree to abbreviate this by writing a = b (mod p). 
Begin by showing that 

Now show that 2S{p)/p = — 3T(p) (modp), where 



np) = E 



We are told that for every integer n with p > n > there is a unique integer m with 
p > m > and m = l/n (mod p). Hence you may now use the result of Problem 1 
to solve Problem 2. 



1.3 Riemann zeta values 

Now let's consider what happens as — )■ oo in the sum S{N) = J2N>n>o ^/n^- The 
answer has a name: ({3), where ({s) is the value of the Riemann zeta function with 
argument s. Bernhard Riemann (1826-1866) found remarkable relations between 
the zeros of ({s) in the complex plane and the distribution of prime numbers. Here, 
we shall study only the Riemann zeta values 

OO -I 

c(^) = E 4 

for integers s > 1. 

Pari-GP is a computer algebra system designed for fast computations in number 
theory and is available from http://pari.math.u-bordeaux.fr for free. It works 
well on Windows and linux platforms. Here, 1 ask Pari-GP to print values of 7c'^/({s) 
for 5 even integers s > 2, at 60-digit precision. 

default (realprecision, 60); 

f orstep(s=2, 10,2, print([s, Pi"s/zeta(s)] ) ) ; 

The output 

[2 , 6 . 00000000000000000000000000000000000000000000000000000000000] 
[4 , 90 . 0000000000000000000000000000000000000000000000000000000000] 
[6 , 945 . 000000000000000000000000000000000000000000000000000000000] 
[8 , 9450 . 00000000000000000000000000000000000000000000000000000000] 
[10 , 93555 . 0000000000000000000000000000000000000000000000000000000] 
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strongly suggests (but does not prove) the exact evaluations 

C(2) 



7r2 



6 

C(4) 



90 

C(6) = 



945 



9450 

C(10) = 



93555 



which were proven in 1735 by the master analyst Leonhard Euler (1707-1783). For 
every even integer s > 2, it is proven that C(s)/vr** is a rational number. 

For folk who love analysis (and only for them) here is a cunning method to prove 
that C(2) = 7r^/6. First one may use binomial expansions to show that 

c(2) = rr^ (3) 



2 Jo Jo 1 + xy 

Adding these results and dividing by 2, we obtain 



Jo 1 — xy 

= rr±^. (4) 



1((2) = /' /' '^^''^ . (5) 

4 Jo Jo 1 — x'^y^ 



Then the delicate transformations 



sm(aj 

X = — (6) 



cos(6) 
sin (6) 
cos(a) 



(7) 



give 



3 7T TT „ 9 

-C(2) = da db 
4 JO Jo 



and hence ({2) = 7r^/6. 



1.3.1 Problem 3: evaluating ^(2) 

Prove equations dS]), (jlj), ([5]) and ([8]). [This is for enthusiasts of analysis, for whom 
hints follow. An apocryphal remark by Sherlock Holmes appears in the solution.] 
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1.3.2 Hints for Problem 3 



4^9 



16 



Show that the integrals in (j3]) and (jlj) give l + | + | + ^ + and 1 
Hence prove (jlj) by considering the difference of these sums. Then work out the 
Jacobian of the transformations and ([7]), namely the determinant 



dx/da dx/db 
dy/da dy/db 



(9) 



and obtain the result in dHl) from 



1.4 Is (^(3)/7r^ an irrational number? 

It is proven that ({3) is not a rational number, but it is not known, for certain, 
whether ({3)/it^ is irrational. It seems, however, overwhelmingly probable that this 
number is irrational. There are rather good algorithms with which to guess integers 
and D such that C(s)/7r* is close to N/D. For oc?c? integers s > 3, no such attempt 
has produced a guess that works. For even integers s > 2, this experimental method 
works rather well. 



1.4.1 Problem 4: the PSLQ and LLL algorithms 

Use the lindep procedure of Pari-GP, at 100-digit precision, to search for an integer 
relation between the three constants in the vector [C(12), tt^^, 1]. Check that the 
PSLQ and LLL options of lindep give the same result and that the resulting formula 
for C(12) still holds at 1000-digit precision. Finally, investigate, with comparable 
vigour, whether there might be an integer relation between the constants in the 
vector [C(13), vr^^, 1]. [This problem is for enthusiasts of free-ware, for whom hints 
follow.] 



1.4.2 Hints for Problem 4 

In the problem, you are asked to perform investigations for C(12) and C(13). Here, 
I indicate how I used Pari-GP to investigate ({2) and C(3)- I began with 

default (realprecision, 100); flag = -3; 
print (lindepC [zeta(2) , Pi^2, 1], flag)~); 

which gave the answer [—6, 1, 0]. This means that, at 100-digit precision, the pro- 
cedure lindep "guessed" the integer relation —6 x ((2) -|-lx7r^-|-0xl = 0. The 
flag —3 invokes the PSLQ algorithm of Helaman Ferguson and David Bailey. Note 
that I introduced the "red-herring" constant 1 as a safety measure. Next, I changed 
the flag to 0, so as to invoke the LLL algorithm of Arjen Lenstra, Hendrik Lenstra 
and Laszlo Lovasz, obtaining the same integer relation. To be even safer, I eval- 
uated 6({2) — vr^ at 1000-digit precision, obtaining the answer O.E-1001 . Then, 
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I investigated the vector [C(3), n^, 1], obtaining no credible integer relation from 
either PSLQ or LLL. The PSLQ algorithm sensibly refused to give me a relation, 
returning instead a large real constant as an error message. The LLL algorithm gave 
me a relation, but this was soon shown to be untenable by increasing the precision. 



2 Sunshine (and Big Bang) numbers 
2.1 The Stefan— Boltzmann constant 

A black body is one that absorbs electromagnetic radiation with maximum efficiency 
and emits it with maximum efficiency, at all frequencies /. It appears black only 
when very cold. When it gets hot it will glow, first red-hot and then at higher 
temperatures T white-hot. The spectrum of its radiation depends on / only via the 
convenient dimensionless ratio 

hf 

where h is Planck's constant and k is Boltzmann's constant. If the temperature, T, of 
the black body (relative to absolute zero) is given in Kelvins (K) and the frequency, 
f, is given in Hcrz (Hz), then the approximate values h — 6.626 x 10"^^ J s and 
k = 1.381 X 10-23 J will serve here. 

The spectrum of sunlight that we receive at the top of the Earth's atmosphere is 
rather similar to the spectrum produced by a black body with T ^ 6000 K. The much 
spikier spectrum that wc receive at sea level results from the complex composition 
of our atmosphere and may be subject to dramatic change, in your lifetime, by our 
carbon emissions. 

Electromagnetic energy is emitted by a black body of surface area A and temperature 
T at rate given by aAT^, where 

^ = ^ ^ = 5.670 X 10-« J m-2 s-^ R-^ (10) 

15 h-^c^ 

is called the Stcfan-Boltzmann constant. It has a dependence on h, k and the speed 
of light c = 2.998 x 10^ m s^^ that is easily found by dimensional analysis. (Just 
look at the units of cr and you can do it in your head.) 

But where, on Earth, does that Stefan-Boltzmann factor 

27r5 

^ 40.80 

15 

come from, for goodness sakes? 

Amusingly, it turns out that we don't need the precise value of this constant to 
relate the temperature of the Earth, Te, to the temperature of the surface of the 
Sun, Ts. If we assume that both the Earth and the Sun radiate like black bodies. 
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then equilibrium, on Earth, is achieved when 



where Re and Rs are the radii of the Earth and Sun and D is the distance between 
them. Note that the first factor on the right hand side is the fraction of the Sun's 
power that reaches us and is determined simply by geometry. Thus we don't need 
the values of a or R^, and arrive at the estimate 



using only the average angular diameter of the Sun, which is about 32 minutes of 
arc. If we take Te = 288 K (about 15 degrees Celsius) as the average temperature 
of the Earth, then we estimate that Tg ^ Te/0.048 = 6000 K. 

2.2 Planck's integral 

In 1900, Max Planck (1858-1947) gave a formula for the black-body spectrum, i.e. 
the probability that the energy emitted by a black body lies in the narrow range of 
frequencies between / and f+df. This is proportional to B{x)dx, with x — hf /{kT) 
and 



When I read about this, in 1962, it was clear that there was an interesting integral 
to do. To derive the factor of 40.80 in the Stefan- Boltzmann constant, one needs to 
prove that 



but the school library only seemed to contain texts saying things like "it can be 
shown that. . . " this evaluation is correct. There was no internet to consult. Nor 
was there a convenient servant like this 

default (realprecision, 60); 

Ihs = intnum(x=0, [[1] ,1] , 2*Pi*x'3/(exp(x)-l) ) ; print(lhs); 
40 . 8026246380375271016988413391247475307067608837433320737990 
rhs = 2*Pi"5/15; print(rhs); 

40 . 8026246380375271016988413391247475307067608837433320737990 

available to a 15-ycar-old in the early 1960's, to lend such numerical reassurance. 

2.2.1 Problem 5: the sunshine number C,{A) 

By expanding in powers of exp(— x) show that 
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2.2.2 Hints for Problem 5 



You will need to prove that 

roo 31 
/ exp(— nx) dx = —. 

J ^ 

To do so, use the integration variable y = nx and integrate by parts, three times. 
2.3 Evaluations of C(4) 

Thus the problem of the power of sunshine has been pushed back to Euler's 1735 
result that C(4) = 7r790. 

2.3.1 C(4) from the cotangent function 

It is a deep result in complex analysis that we have the wonderful formula 

^= £ — . (11) 

The right hand side has the same singularities, aX z = nn, with the same unit residues 
as occur for the cotangent, on the left hand side. Thus the difference between the 
right and left hand sides is an entire function, with no singularities in the complex 
plane. That does not, of itself, mean that this difference is the zero function, since 
exp(2;) is an example of an entire function that does not vanish. Let's take ( ITTj) on 
trust, from deeper thinkers, and use it to show that ({4) = 7r^/90. First we should 
multiply by z, to remove the singularity at 2 = 0, and then combine the terms with 
positive and negative n, to obtain 

^cosf^) _ o 2 1 



sm[z ) n?'n'^ — z^ 



with a sum that converges near 2; = 0. We know how to expand the left hand side, as 
z — )■ 0, using a ratio of very simple Taylor series for cos(z) and sin(2;). On the right 
hand side we find the constants C,{2k) / ti'^^ , after binomial expansion. In particular, 
by simplifying the relation 



z'^l2\ + z^/A\ + 0{z^) . z- z- 



v2 ^4 



l-2C(2)--2C(4)- + 0(^' 



I- z^/V. + z^/h\ + 0{z^) ^''n^ 'tt 

we easily prove Euler's results that ({2) = and C(4) = vr^/QO. We may go 

further, if we have a convenient servant to do the Taylor expansions for us: 

print(Vec(l - z*cos(z)/sin(z) + 0(z"14))/2); 

[1/6, 0, 1/90, 0, 1/945, 0, 1/9450, 0, 1/93555, 0, 691/638512875, 0] 
shows the origin of Euler's prime 691 in the numerator of ({12) /tt^"^. 
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2.3.2 ^(4) from Fourier analysis 



Here I shall be even sketchier, since the context is more strongly related to my 
undergraduate studies than to my school work. Joseph Fourier (1768-1830) provided 
a method for representing a function, say /(x), on a finite interval, say it > x > — tt, 
by an infinite series of sinusoidal functions, say cos(?t,x) and sin(na;), with coefficients 
that may be evaluated from integrals of the product of f{x) and the sinusoidal 
functions. His method looks much prettier if we use Euler's famous formula 

exp(iz) = cos(z) + isin(z) 

where i^ = —1. But don't be put off by "the square root of —1"; we shall not need 
it, to evaluate C(4)- 

There are not many questions on Fourier analysis that can be solved straightfor- 
wardly, under examination conditions. Normally, the student is asked to analyze a 
very simple function, say f{x) = x^, for which the Fourier coefficients in the series 

oo 

fix)= H anexp(inx) (12) 

n=— oo 

may be computed fairly easily. Now suppose that the student obtains the correct 
Fourier coefficients, namely 

_ j-^ exp{-mx) X , _ r cos{nx) 2 i _ / ^r^/S forn = 

""""l. 2^ }{x)dx-J^ dx - j 2(_i)n/^2 forn^O 

then s/he might be asked to substitute x = in the Fourier series ( |T2l) and obtain, 
thereby, an evaluation of the Riemann eta value 77(2) = Z]n>o(~l)"~^/''^^ = C(2)/2. 
This too is rather easy: at x = we obtain = vr^/S — 4?7(2) from the Fourier 
series (fT2|) . using the result for the Fourier coefficients. 

An interesting exam question might say: "use Parseval's theorem to evaluate C(4)"- 

The theorem in question gives an integral of the square of any function /(x) as the 
sum of the squares of its Fourier coefficients a„: 



- / \f{x)\'dx= 



2.3.3 Problem 6: C(4) from Fourier and Parseval 

Use the given Fourier series for x^ and Parseval's theorem to prove that ({4) = 7r^/90. 
2.4 The sunshine number ({3) 

We have seen that a black body of surface area A and absolute temperature T 
radiates energy at a rate 

energy ^ {kT)*A 
time ^ h^c^ 
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where 

derives from an integral over the black-body spectrum. 

Like all forms of electromagnetic radiation, hght is emitted and absorbed as "bun- 
dles" of energy called photons. The rate at which a black body of surface area A 
and temperature T emits photons is 

photons _ ^ {kTfA 
time ^ h^c^ 

where 

^3 = / r\—{ = 47rC 3 ^ 15.11 

Jo exp(x) — 1 

2.4.1 How many photons in your oven? 



volume y tic J 

where the constant is now 4/3 and the factor of 4 takes account of the fact that only 
\ of the photons are moving towards the nearest wall and these have an average 
value of I for the cosine of the angle of incidence with which they hit the wall. 

In the U.K. the formula for the absolute temperature of an oven at "gas mark" 
M > 1 is 

T{M) = (^273.15 + - 1) + 135^ K 

since M = 1 is defined as 135 degrees Celsius (identical to 275 degrees Fahrenheit) 
and then each new mark increases the temperature by 25 degrees Fahrenheit, i.e. by 
125/9 degrees Celsius. The zero of the Celsius scale is defined as 273.15 K, which is 
very close to the melting point of pure ice at standard atmospheric pressure. 

Let's set an empty 40-litrc oven at gas mark 9 and assume that the walls and door are 
perfectly black, on the inside. After this oven has arrived at the design temperature, 
we switch off the power. Working out 

A;T(9) 1.381 x lO^^s x 519.3 1 . , 4 1 

— ^ = m-^ = 3.61 X 10^ m~^ 

he 6.626 X 10-34 X 2.998 x lO^ 

we expect to have 

167rC(3) X (3.61 x 10^)^ x 0.040 = 1.14 x 10^^ photons 

inside the oven. When the oven eventually cools down to room temperature, with 
T 300 K, the number will have decreased by a factor of (300/519)^ = 0.193, so 
now there will be merely 22 trillion photons inside. 

Now let's turn down the oven to 2.73 K. 
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2.5 How many photons in your Universe? 



About 13.7 billion years ago, the whole of the Universe was at a temperature com- 
parable to the present temperature of the surface of the Sun. Photons, electrons, 
protons (and other particles) were in constant collision. But then the expansion of 
the Universe (in fact the expansion of space itself) cooled things down and when 
the temperature fell below about 3000 K the photons ceased to interact significantly 
with matter. Since then, space has expanded by a factor of about 1000, in each of its 
three dimensions, and the cosmic background radiation that we readily detect nowa- 
days has a black-body spectrum with a temperature of about 3 K. The expansion of 
space has stretched out the most probable wavelength, A = (2.898 x 10~^ m K)/T, 
taking it from the visible region, at the temperature of the surface of the Sun, to the 
microwave region, now. In the last few years we have come to understand a great 
deal about the hot Big Bang in which our Universe originated, by studying this 
cosmic microwave background radiation. Its present temperature is known rather 
precisely: 

To = (2.728 ±0.002) K 
where the subscript means "now" . 

The whole of the Universe now has a background photon density 

cosmic photons ( kTn\^ , _ „o , , „, 

— f— = 167rC3 — ^ =4.12xl0'm-l 13 

volume of Universe \ nc J 

What "volume" should we multiply by, to count the total number of photons? For 
all we know, the Universe might have an infinite volume. Yet we cannot see all of 
it. Let's consider just that volume 

out to the present distance Rq of the places from which we now receive light from 
the Big Bang. Further than that we cannot see. It is believed, nowadays, that wc 
are allowed to use Euclidean geometry in this calculation. Space-time is curved but 

space seems not be, on a cosmic scale. 

You might think that Ro/c ought to be close to to = 13.7 x 10^ years, the present 
age of the Universe. In fact, it's not as simple as that: the curvature of space-time 
in the early Universe means that we must use the general theory of relativity. In 
the simplest current model, 

-Rq _ A(^A,o) 
cto hi^hfi) 

is the ratio of two integrals that depend on the crucial number Vt\fi = 0.734 ±0.02, 
which is the fraction of the present space-time curvature of the Universe that comes 
from the cosmological constant. (You may have heard this contribution referred to, 
rather unhelpfully, as "dark energy".) The integrals are 



^a^c+{l - c)a 
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^0 



a da 1 / 1 + ^/c^ 

m 



and the first is too hard to do by hand. So let's ask our servant to do it: 

default (realprecision, 5) ; c = 0.734; s = sqrt(c); 
print (intnuiii(a=0,l, l/sqrt(a'4*c+(l-c)*a))) ; 
3.4700 

print(log((l+s)/(l-s))/3/s); 
0.99677 

There are about are 3.156 x 10^ seconds in a year. So hght travels 3.156 x 10^ x 
2.998 X 10^ = 9.46 x 10^^ metres in a year. Hence I claim that there are about 

47r /_3^ ^ ^ Q ^ ^ ^ i5\ ^ X 10^ = 1.58 x 10^^ photons 
3 V0.997 J ^ 

present within the horizon beyond which you cannot possibly look. 

Please remember that we needed the Riemann zeta value 

OO 1 

C(3) = ^ — = 1.2020569... 



n=l ^ 



to do this cosmic calculation. 



2.6 How many neutrinos in your Universe? 

In nuclear beta-decay a neutron turns into a proton. An electron and a anti-neutrino 
are ejected from the nucleus: 

n — > p + e~ +Ve 

sharing the released energy. We call Ve the "electron anti- neutrino" . Other nuclei 
eject positrons and electron neutrinos: 

p — > n + e'^ + 1^^ 

A second pair of neutrinos is produced by the decays 

A*" — > e~ i/^ 

jj.'^ — > e+ -M7^ i/e 

which occur in the Earth's atmosphere, where muons are created by cosmic rays. A 
third pair, Vj- and i^^, is readily produce by the decays of "tau-leptons" , r^, at high 
energy physics laboratories. There seem to be no more. 

There is a cosmic background of neutrinos from the hot Big Bang with a ratio 

neutrinos ^ Jg j^^^ ^-^^-^ 



photons h \ Tq 
12 



of the neutrino and photon densities, where ^ 1.95 K is the present temperature 
of the neutrino background and 



with a smaller integral J3 for the neutrinos, because they are particles that obey 
an "exclusion principle" . Neutrinos, like electrons, utterly refuse to share quantum 
states with particles identical to themselves. (This rather selfish behaviour is fortu- 
nate for us, because it is the exclusion principle for electrons that is responsible for 
the wonderful chemistry revealed by the periodic table of the elements.) 

The factor of 3 in the ratio f lT^ comes from the three "flavours" of neutrino: e, fi 
and r. The next factor comes from the exclusion principle. In the simplest model, 
the final factor is a rational number: {Tjj/TqY = ^• 

2.6.1 Problem 7: C(3) integral for cosmic neutrinos 

Show that J3 = 37rC(3) and hence that there are, on average, about 337 neutrinos, 
left over from the Big Bang, in each cubic centimetre of your Universe. [Neutrinos 
permeate matter with great ease. There may be more than 300 cosmic neutrinos in 
each cubic centimetre of your brain, as you solve this problem, along with higher 
energy neutrinos from local sources, such as the Sun and decays of muons in the 
Earth's atmosphere.] 

3 Feynman's zeta values 
3.1 Particles and fields 

On January 25, 1947, the journal Nature published a letter from Donald Perkins [13] 
(ne 1925), entitled "Nuclear disintegration by meson capture". 

On December 30, 1947, the Physical Review received a one-page article by Julian 
Schwinger [15] (1918-1994) entitled "On quantum-electrodynamics and the magnetic 
moment of the electron" . 

Between January and December of that year came momentous discoveries of two 
new types of particle, muons and pions, and, quite independently, the stimulus for 
great leaps of understanding in the quantum field theory of electromagnetism, by 
Schwinger, Hans Bethe (1906-2005), Richard Feynman (1918-1988) and Freeman 
Dyson (ne 1923). 
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3.1.1 Particles 

Pions (or pi-mesons as they were called in 1947) are created when high-energy pro- 
tons from distant regions of the Universe hit the Earth's atmosphere. A negatively 
charged pion decays to give (most often) a muon and a neutrino: 

The muon then decays to give an electron and a pair of neutrinos. The discovery of 
pions and muons, in 1947, by the tracks that they leave in photographic emulsion, 
was the beginning of the experimental side of my subject: particle physics. Since 
then many more types of particles have been found, in collisions at particle acceler- 
ators. Currently, experimenters at the large hadron collider in Geneva are searching 
for a theoretically expected particle, called the "Higgs boson" . 



3.1.2 Fields 

James Clerk Maxwell (1831-1879) described the electromagnetic fields, E and B, 
that are created by, and in turn influence, charged particles. An example is the force 
between a pair of current carrying wires, which is used to define the Ampere (A), as 
a unit of electrical current, in terms of the definition of the Newton (N), as a unit 
of force. This definition appears in Maxwell's equations via a constant 

^0 = i!L N A-2 

which appears in magnetic-field calculations. Using this value, we may define a very 
important dimensionless number 

e^c 1 



2h 137.0359990 . . . 

called the fine structure constant. It is know to an accuracy of about one part in a 
billion. Here is the square of the electron's charge, c is the speed of light and h 
is Planck's constant. Physicists tend to remember the reciprocal of the value of a, 
which is close to 137. 

The quantum field theory developed in the late 1940's, by Schwinger, Feynman 
and Dyson, combined Maxwell's classical theory of clectromagnetism (in which 
occurs) with Einstein's special theory of relativity (in which c is rather important) 
and earlier quantum theory (in which h is rather important). 

Werner Heisenberg (1901-1976), Erwin Schrodingcr (1887-1961) and Paul Dirac 
(1902-1984) had developed the theory of quantum mechanics, between 1925 and 
1932. This confronts the fundamental unpredictability of individual events, at the 
deepest level of physics, and gives ways of calculating the probabilities of outcomes, 
in a large number of measurements. But it does not handle the creation and annihi- 
lation of particles. Schwinger and Feynman solved that problem, in the interactions 
of photons, electrons and positrons. 
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In the 1930 's it was realized that a "photon field" was needed, to explain the creation 
of photons by the interactions of charged particles, and that an "electron field" was 
needed, to explain the creation of electrons and positrons by the photon field. 

I shall try to indicate how the "sunshine" number C(3), and its cousins C(5), CC''), 
etc, appear in the application of quantum field theory to particle physics. But first, 
it would be a good idea to look at the result of a rather simpler calculation. 



3.1.3 Electron-positron creation in pion decay 

The positron, e^, is the anti-particle of the electron, e~, with exactly the same mass, 
me, and exactly the opposite charge, e. 

The neutral pion 7r° usually decays into a pair of photons: 

7r° — > 7 + 7 



where 7 is the symbol that we use for a photon (sometimes called a gamma ray). In 
the rest frame of 1 
is the mass of vr''. 



the rest frame of the pion, each photon carries away an energy |m^oc^, where m^o 



In 1951, Richard Dalitz (1925-2006) used the quantum field theory of Feynman and 
Schwinger to predict another, rarer decay: 

7r° — ^ 7 + e" + e+ 

in which an electron-positron pair is created, in place of one of the photons. His 
predicted rate was a fraction J9j 

Dalitz pair probability = — (-In ( — ~'~] — 1-185% 



TT V3 V rric J 3 

of the total decays. This was a rather successful prediction. The best result from 
modern experiments is that (1.198 ± 0.032)% of the decays of vr" have the electron- 
positron pair predicted by Dalitz. 

Dalitz had to do an integral to make his prediction, since the total energy released 
is here shared between 3 particles, in the final state, and we must integrate over all 
the possible ways of doing that. In any individual decay, we cannot say how the 
energy will be shared. We can predict only the probabilities of the various outcomes 
and check those predictions against a large number of measurements. 

3.2 Magnetic moment of the electron 



magnetic moment , a , f ct\^ , f Ci\^ 

I + A1- + A2 - +A3 - + 



Bohr magneton vr Vtt/ Vtt 
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where the Bohr magneton 



- 9.274 X 10"^^ J T"^ 



Annie 

has the units of Joules per Tesla and is the value predicted by Dirac in 1928. In 
1947, Schwinger found the first correction term, with [15] 

4, = i 

^ 2 

In 1950, Robert Karplus and Norman KroU claimed the value [10] 

125 2687 
28C(3)-54C(2)ln(2) + — C(2)-^ = -2.972604271... 

for the coefficient of the next correction. It turned out that they had made a mistake 
in this rather difficult calculation. The correct result [161 [E 
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A2 = -C(3)-3C(2)ln(2) + -C(2) + — = -0.3284789655... 

was not obtained until 1959. Not until 1996 was the next coefficient found 

As = -^C(5) + ^C(3)C(2) - fc(4) - f f/3,1 

139,, , 596,, , , 17101,, , 28259 , , 

+ C3 C2 n2 + ^2 + 15 

= 1.181241456... 

by Stefano Laporta and Ettore Remiddi. The irrational numbers appearing on the 
second line are those already seen in A2. On the first line we see zeta values and a 
new number, namely the alternating double sum 

U3i= y ^ { ^ -0.1178759996505093268410139508341376187152 . . . 

I visited Stefano and Ettore in Bologna when they were working on this formidable 
calculation and recommended to them a method of integration by parts, in D dimen- 
sions, that I had found useful for related calculations in the quantum field theory of 
electrons and photons. Here D is eventually set to 4, the number of dimensions of 
space-time. But it turns out to be easier if we keep it as a variable until the final 
stage of the calculation. Then if we find parts of the result with factors like 1 / (D—A), 
or 1/{D — 4)^, we need not worry that these parts tend to infinity as D — > 4; all 
that matters is that the complete result is finite. Based on my D-dimensional expe- 
rience [6j, I expected their final result to look simplest when written in terms of the 
new number U^^i. This indeed turned out to be case. If, however, one writes (|T5|) in 
terms of a less appropriate non-zeta value, such as [5] 

°° 1 / 1 \ " 1 1 1 1 

£ (2) = -24l'"P" + + - 

then the additional irrational numbers [ln(2)]^ and C(2)[ln(2)]^ will also appear, since 
they I had expected, absent when one uses f/3^1. 
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3.3 Feynman diagrams 

X 
15 




Every Feynman diagram tells a possible story. In this one, a muon emits photon 
1, which creates an electron and a positron. These exchange photons 2 and 3 and 
annihilate to give photon 4. Meanwhile, the muon interacts with a magnet X, via a 
(very low energy) photon 5, and then the muon neatly absorbs photon 4, to emerge 
from this story with (almost) the same energy and momentum as it started with. 

This Feynman diagram specifies only one of the very large number of integrals that 
one must do to compute the term of order in the magnetic moment of the muon. 
Feynman found methods that give this particular contribution as an 8-dimensional 
integral. 

In 1995, Pavel Baikov and I found a method of obtaining a numerical accuracy of 
5 significant figures for such integrals pLj. A previous method, from Toichiro (Tom) 
Kinoshita [11], was shown to give errors of 5%. 

Tom is one of my true heros. He was born in 1925 and studied with Sin-Itiro 
Tomonaga (1906-1979), who shared the 1965 Nobel Prize in physics with Schwinger 
and Feynman. 

3.4 Multiple zeta values 

When there is a single very large external energy scale, it is often a good approx- 
imation to neglect the masses of particles inside a Feynman diagram. This occurs 
in the case of electron-positron collisions at very high energy, which produce large 
numbers of particles, such as protons and pions, containing light quarks. So far 
all the results obtained from such diagrams have been zeta- valued, with (^(7) now 
appearing in recent results [2]. 

I believe that this cannot continue, indefinitely. In 1995, Dirk Kreimer and I calcu- 
lated some diagrams in a "toy" theory that has only one type of interacting quantum 
field, unlike the more pertinent electron-positron-photon and quark-antiquark-gluon 
theories of the electromagnetic and strong interactions. In this toy theory, we found 
Feynman diagrams that give a multiple zeta value [HI [7J H] 

C(5,3)= E ^ 
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which is the first double sum that cannot be reduced to zeta values. 



So I end with the conjecture that when practical quantum field theory progresses 
just one more step in difficulty, to calculations that yield the zeta value C(9); then 
those same calculations will also yield the multiple zeta value C(5, 3). 



4 Solutions 
4.1 Solution 1 

Let f{x) = -^x{x — l){2x — l)(3a;^ — 3a; — 1). Then simple algebra shows that 
f{x + 1) — f{x) = x^. Hence it follows, by induction, that Z]Ar>n>o'^^ = f{^): 
each integer N >2, since this is true for N = 2, 



4.2 Solution 2 

We replace n by j9 — n in the summation S{p) = J2p>n>o ^/n^, add the result to S{p) 
and divide by the prime p to obtain 2S{p)/p = J2p>n>o s{p, n), where 

, , 1/1 1 \ p^-3pn + 3n^ 3 , ^ , 

s p, n = - — + = -— = mod p). 

p \rf^ [p — n)-^ J n-^\j) — ny^ 

Hence 2S{p) = (mod p) and 2S{p)/p = —3 J2p>n>o (mod p). For each integer 
n with p > n > there is a unique integer m with p > m > and m = 1/n (mod p). 
Hence we replace n by m and obtain 2S{p)/p = —3 J2p>ni>o (mod p). Finally we 
use the result ([2]) of Problem 1, which shows that 30 X]p>m>o "^'^ = (mod p) and 
hence that 20S{p)/p = (modp). Thus the numerator of S{p) is divisible by for 
every odd prime p that is not equal to 5. 



4.3 Solution 3 

For \z\ < 1 we may use the binomial expansion 1/(1 — 2;) = J2n>o^"'~^- With 
z = xy under the integral sign in (|3]), this gives J2n>o ^/n"^ = C(2), as claimed. With 
z = —xy in it gives J2n>oi~^)"'~^ / ■ C!all this latter sum 77(2) (a Riemann 
eta value). Then ((2) - r/(2) = 2Efc>ol/(2fc)^ = |C(2) and hence r/(2) = ^((2), as 
claimed in (jl]). Now we compute the Jacobian in Qj, which gives 

^ cos(a) cos(6) sin(a) sin(6) sin(a) sin(6) 
' cos(6) cos(a) cos^(a) cos^(6) 

= 1 — tan^(a) tan^(6) = 1 — x^y^. 

Hence we may replace dxdy/{l — x'^y'^) in ([5]) by the infinitesimal product dad6. 
Finally we have to take care of the limits for a and h. We must impose the condition 
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a + b < ^ to ensure that neither x nor y exceeds 1. Hence we obtain half of the area 
of a square with side |, in ([8]), which proves that ({2) = 

[This is called an "elementary" proof. As you have seen, this does not mean that 
it is an "easy" proof. Rather, it means that the proof does not rely on some deep 
result in complex analysis. Later, you will see proofs that may appear to be much 
easier, yet rely on much deeper assumptions. When Sherlock Holmes (in the movies, 
if not in Conan Doyle's books) declares that his line of reasoning is "elementary, 
my dear Watson" , he does not mean that it is trivial; rather he means that he can 
explain it, sometimes at length, without relying on external authority. The same 
type of thing occurs in mathematics. Frits Beukers, Eugenio Calabi and Johan Kolk 
published this "elementary" proof in 1993 |3].] 



4.4 Solution 4 

For C(12), both PSLQ and LLL should have given you Euler's result 



' 638512875 

with the intriguing prime 691 in the numerator. You should not, according to 
current belief, have found any significant integer relation between the constants 
[C(13), vr^^, 1]. If you believe that you did, then please check it at higher precision, 
where it will almost certainly evaporate. 



4.5 Solution 5 

The black-body function may be written as B{x) = 27ra;^ Z]n>o 6xp(— nx). To evalu- 
ate I4 = /o°° B{x) dx we make the substitution x = y /n, in the n-th term, obtaining 
/4 = 27rr(4)(^(4), with the term in ^(4) = J2n>o 1/"^^ coming from the change 
of variables and an overall factor of 

roc 

r(2;) = / i/^"^exp(-2/) dy 
Jo 

appearing at z = 4. Integration by parts gives T{z) = {z — l)T{z — 1). Moreover, 
r(l) = 1. Hence, for positive integer n, we obtain T{n) = {n — 1)! after n — 1 
integrations by parts. Thus I4 = 1271(^(4). 



4.6 Solution 6 



With f{x) = x^, Parseval's theorem gives 7i'^/5 = (7r^/3)^ + 8^(4) and hence ({4) = 
(1/5 - l/9)7r78 = 7r790. 
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4.7 Solution 7 

The expansion of l/(exp(x) + 1) — X]n>o( — 1)" ^exp( — nx) gives a Riemann eta 
value in J3 = Anrj^S). Subtracting 77(8) = J2n>oi~^)"'~^ / ''^^ from ^(3) = J2n>o l/""-^, 
we obtain 

C(3)-.(3)^2f pLj4c(3). 

Thus J3//3 = r7(3)/C(3) = | and the neutrino to photon ratio is 3 x | x ^ = ^. We 
know from ( !T2|) that there are about 412 photons per cubic centimetre. Multiplying 
by we obtain 337 neutrinos per c.c. 



5 Acknowledgments 

5.1 The Manchester Grammar School: 1958—1964 

I learnt calculus from Eric Hodge (1904-1962) who sadly died at the end of my 
fourth-form year. My enthusiasm for the method of induction comes from him. The 
somewhat awesome High Master of MGS at that time was Eric James (1909-1992) 
who had taught Freeman Dyson at Winchester. Dyson's latest book, The Scientist 
as Rebel (2006), is dedicated to Eric and Cordelia James. 

My attention to detail in mathematical calculations was nurtured in the sixth form 
by Philip Scofield (1930-2010) for whom a sharpened piece of chalk was the most 
precise tool in all analysis. My physics teacher Aden Womersley (1931-2007) had 
a lively interest in modern developments and knew about quarks as soon as Gell- 
Mann hypothesized them in 1964. Yet my most influential teacher was John Scobell 
Armstrong (1927-2001) whose love of enlightenment philosophy spoke to both my 
heart and my head. 

Around the time of the Cuban missile crisis of late 1962, as I recall, I became 
fascinated by the large factor 27r^/15 ~ 40.80 in the Stefan-Boltzmann constant that 
determines how much sunshine hits the upper atmosphere of planet Earth. Trying to 
understand its origin seemed like a good diversion from wondering whether Kennedy 
and Khrushchev were about to release catastrophic amounts of nuclear energy. 



5.2 Zambia: 1965 

At the end of 1964, I left the Manchester Grammar School and taught for 9 months 
at Munali Secondary School, in Lusaka. The discipline of teaching physics and 
maths to bright and highly motivated students did me a power of good. One of 
my former students, the late Dr Wedson C. Mwambazi, became the World Health 
Organization representative for Tanzania, and later for Ethiopia. 
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5.3 Oxford: 1965-1968 

I was very fortunate to have Heinrich Kuhn (1904-1994), David Brink (ne 1930) 
and Patrick Sandars (ne 1935) as my physics tutors at Balhol Cohege. From Don 
Perkins' lectures I discovered particle physics and became determined to study this 
subject as a post-graduate student at the newly established University of Sussex. I 
also heard in Oxford how seriously Dick Dalitz regarded the quark model. 

5.4 Sussex: 1968-1971 

I was equally fortunate in my postgraduate teachers, Gabriel Barton, David Bailin 

and Norman Dombey. It was fun, in my first term, to work out how Schwinger 
had obtained the correction term i(a/7r) for the magnetic moment of the electron, 
back in 1947, the year of my birth. Stan Brodsky visited Sussex and told me how 
amazingly difficult it would be to obtain an exact value for the coefficient of (a/vr)^. 

5.5 Stanford: 1971-1973 

This good fortune continued during my post-doctoral years. The Stanford Linear 

Accelerator Center (SLAC) in California proved just the right place to be, when 
physicists were discovering the "parts" of the proton, by scattering electrons off 
protons. Sid Drell and Jim Bjorken were great enthusiasts of Feynman's "parton 
model". Ken Johnson (1931-1999) was visiting SLAC. He was already sure that 
quantum field theory had come back to stay and would soon make sense of the 
strong interactions between pions and protons. 

Yet it was also interesting to hear of competing metaphysics, from Geoff Chew at 
Berkeley, across the Bay. Geoff thought that there might be no quantum field theory 
of strong interactions. Instead, he advocated a "bootstrap" approach that might free 
us from previous reliance on the "undemocratic" supposition of forces between a few 
types of fundamental particle, somewhat in the manner of Baron von Miinchhausen, 
who by legend tried to defy gravity by pulling himself up by his own boot laces. 

At Stanford, however, we desperately wanted Feynman's partons to have the frac- 
tional charges and the "three mathematical colours" that had notionally been as- 
signed to quarks. Those were the glory days, when experiment was the ruler of our 
subject, and we had to wait but a short while for extra data, from electron-positron 
collisions and neutrino-proton collisions, to confirm the quark-parton model and 
open the fiood gates to a new quantum field theory. 

5.6 Geneva: 1973-1974 

At the laboratory of the European Council for Nuclear Reseach, CERN, in Geneva, 
I learnt from Gerard 't Hooft how quantum field theory could also make sense of the 
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weak interaction, responsible for muon decay. Gerhard helped me to understand 
why strong interactions get weaker and weak interactions get stronger, at higher 
energies. Ken Johnson was right: quantum field theory is definitely here to stay. 

I came to appreciate just how much hard work there was to do to meet the computa- 
tional challenges of computing the Feynman diagrams of these new theories. It was 
Gerhard's supervisor, Martinus (Tiny) Veltman, who had begun to confront these 
challenges in the early 1970's. In 1974 there appeared a review of Tiny's computer- 
algebra package, Schoonship, whose modern-day successor Form has been developed 
by Jos Vermaseren to be a powerful engine which implements the subtle algorithms 
that yield zeta-values like C(7) in current calculations of massless Feynman diagrams. 

5.7 Oxford: 1974-1975 

Returning to Oxford as a junior research fellow, I revelled, with Dick Dalitz, as 
wonderful new results from electron-positron annihilation came from Stanford, and 
other laboratories, making the case for quarks more and more compelling. 

One day, when I was scheduled to give a colloquium on these new results, a school- 
boy from Eton came to grill me. His grasp of modern physics was remarkable 
and I wondered what he might later contribute to my subject. In fact, after a 
brief but distinguished period of work, Stephen Wolfram forsook particle physics. 
His programme Mathematica has been marketed with great commercial vigour and 
is now widely used by students and researchers to perform some of the simpler 
computations in particle physics. 

5.8 The Open University: 1975— present 

When I joined the Open University, one of my first tasks was to study Steven 
Weinberg's fine book. Gravitation and Cosmology, and then help to write a course, 
Understanding Space and Time. Julian Schwinger came to Milton Keynes, to work 
with us on producing this course. At my earnest request, he gave us a memorable 
3- hour extempore account of general relativity, as he perceived it. I could scarcely 
believe my luck. I was getting paid to learn things, for myself, and then was able 
to help to share those good things with open-minded and enthusiastic students. My 
own openness of mind is reinforced by my avocation, by my partner Margaret and 
by our sons Stephen and Peter. 

5.9 Envoi 

To all those acknowledged above, and to many more, I give humble and hearty 
thanks for all their goodness and loving kindness. 
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